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ABSTRACT 

We present observations made with the Berkeley- Illinois-Maryland Association millimeter array of the 
H2S 2(2,0)^ 2(1,1) and C^*0 2^ 1 transitions toward a sample of four hot molecular cores associated 
with ultracompact Hii regions: G9.62+0.19, GlO.47-1-0.03, G29.96-0.02 and G31.41-h0.31. The angular 
resolution varies from 1.5 to 2.4 arcsec, corresponding to scales of ~0.06 pc at the distance of these 
sources. High-velocity wings characteristic of molecular outflows are detected toward all four sources 
in the H2S line. In two cases (G29.96 and G31.41) red- and blueshifted lobes are clearly defined and 
spatially separate, indicating that the flows are collimated. We also confirm the previous detection of the 
outflow in G9.62F. Although the gas-phase H2S abundance is not well constrained, assuming a value of 
10~^ yields lower limits to total outflow masses of ~8 Mq, values which are consistent with the driving 
sources being massive protostars. Linear velocity gradients are detected in both C^*0 and H2S across 
G9.62, G29.96 and, to a lesser extent, G31.41. These gradients are observed to be at a different position 
angle to the outflow in G9.62F and G29.96, suggestive of a rotation signature in these two hot cores. Our 
observations show that these hot cores contain embedded massive protostellar objects which are driving 
bipolar outflows. Furthermore, the lack of strong centimeter-wave emission toward the outflow centers 
in G29.96 and G31.41 indicates that the outflow phase begins prior to the formation of a detectable 
ultracompact Hii region. 

Subject headings: stars: formation — HII regions — ISM: molecules — radio lines: ISM 



1. INTRODUCTION 

Hot molecular cores are compact, dense cores of gas 
which display a rich array of chemical species and are fre- 
quently found in association with ultracompact Hii re- 
gions (e.g. Kurtz et al. 2000). Temperatures in the re- 
gion of 100-300 K have been derived from ammonia and 
methyl cyanide observations (Cesaroni et al. 1994, Olmi et 
al. 1996). In addition to the high temperature a striking 
feature of hot cores is the elevated abundance of saturated 
molecules and rare species compared with dark clouds such 
as TMC-1, which chemical models show can only arise if 
chemical reactions on grain surfaces are included (Brown, 
Millar & Charnley 1988). However, it has been unclear as 
to whether all hot cores marked the actual sites of massive 
star formation or if some were remnant clumps heated by 
nearby OB stars (e.g. Watt & Mundy 1999). 

Many studies have been made of hot cores, and observa- 
tions with interferometers have begun to reveal the geome- 
try of the gas and dust in these regions (Blake et al. 1996; 
Wright, Plambeck & Wilner 1996). Small-scale velocity 
gradients have been inferred from methyl cyanide obser- 
vations (Olmi et al. 1996) but thus far the interpretation 
is ambiguous. Subarcsecond observations of the (4,4) in- 
version transition of ammonia by Cesaroni et al. (1998, 
hereafter C98) represent the highest resolution molecular- 



line images published to date, resolving the hot cores in 
G10.47-t-0.03, G29.96-0.02 and G31.41-H0.31. C98 also 
observed temperature and velocity gradients in these three 
sources and concluded that the hot cores probably encom- 
pass rotating flattened structures, perhaps disks. 

Further evidence for the presence of embedded sources 
within hot cores came from Testi et al. (2000) who discov- 
ered a faint 3.5-cm radio continuum source within G9.62F, 
while Hofner, Wiesemeyer & Henning (2001) detected 
HCO^ line wing emission with a spatially bipolar distri- 
bution, interpreted as a molecular outflow, in the same hot 
core. Beltran et al. (2004) detected a velocity gradient in 
G31. 41-1-0. 31 which they interpreted as rotation about a 
central source driving a bipolar flow (seen in ^^CO: Olmi 
et al. 1996). 

In order to shed further light on the nature of hot 
cores, we have begun a high-resolution study with the 
Berkeley-IUinois-Maryland Association (BIMA) millime- 
ter array of a sample of four hot core sources: G9.62-f 0.19, 
G10.47-f0.03, G29.96-0.02 and G31.41-h0.31 (hereafter 
G9.62, G10.47, G29.96 and G31.41 respectively). In this 
paper we present observations of the H2S 2(2,0)^2(1,1) 
and C^^O 2^1 lines in order to explore the kinematics 
of the gas associated with the hot cores. In a companion 
paper, Wyrowski, Gibb & Mundy (2004, hereafter Paper 
II) present observations and modeling of the continuum 
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emission at 2.7 and 1.4 mm with sub-arcsecond resolu- 
tion of these objects. They find that in all cases but 
G10.47 the millimeter emission is clearly offset from the 
free- free emission seen at centimeter wavelengths and that 
the strong 1.4- mm continuum emission is evidence that 
these hot cores are massive, dense cores. 

H2S is one such saturated species mentioned above, only 
forming on dust grains or by hydrogenation reactions in 
hot shocked gas (e.g. Charnley 1997). Models of the gas- 
phase chemistry involving H2S predict that over a period 
of ~10'' yr H2S reacts with oxygen-bearing species to form 
the daughter products SO and SO2, causing the gas-phase 
abundance of H2S to decline from ^ 10~^ to 10""'^° or 
lower (e.g. Rodgers & Charnley 2003; Wakelam et al. 
2004). Thus its distribution will trace the regions where 
dust grains have been recently subject to sufficient heat- 
ing to evaporate material from the grain mantles (Hatchell, 
Roberts & Millar 1999), cither thermally by the radiation 
from nearby OB stars or non-thermally through the action 
of shocks in an outflow. The gas phase abundance of H2S 
is not well known, with a number of studies giving values 
in the range ~ 10"^° to 10"'^ (van der Tak et al. 2003; 
Minh et al. 1991). H2S has been observed in the outflow 
associated with the Orion hot core by Minh et al. (1990) 
who derive abundances as high as a few times 10~° rela- 
tive to H2. As we show below, we detect high- velocity H2S 
emission which also originates in outflows from sources em- 
bedded within the hot cores. 

2. OBSERVATIONS 

The BIMA array (Welch et al. 1996) was used in A (ex- 
cept G9.62), B, C and D conflgurations on various dates 
between 1999 May and 2001 February. Nine antennas of 
the array were employed, tuned to a lower-sideband fre- 
quency of 216.7 GHz (corresponding to a wavelength of 
1.4 mm). The intermediate frequency of 1.5 GHz places 
the upper sideband at 219.7 GHz. The correlator setup 
included one spectral window tuned to the frequency of 
the H2S 2(2,0)^2(1,1) at 216.7104 GHz with C^^O 2^1 
at 219.5603 GHz falling in the upper sideband. The ob- 
serving bandwidth was 50 MHz, with a spectral resolution 
of 781 kHz (corresponding to 68.7 and 1.07 kms~^ respec- 
tively). 

Antenna gains were checked with observations of plan- 
ets and found to be accurate within 20 per cent. Phase 
calibration involved observation of the quasars 1733—130, 
1743-038 and 1751-^096 for 3-4 minutes every 15-30 min- 
utes (depending on array configuration). Observations 
made in the A configuration (as well as the B configu- 
ration for G9.62 and G29.96) also utilized fast switching 
between source and calibrator every 1.5-2.5 minutes to 
track the phases (see Looney, Mundy, & Welch 2000 for 
further details). The continuum emission was subtracted 
from the visibility data using channels free from line emis- 
sion. The contimmm emission from each source at 1.4 mm 
is strong enough to permit self-calibration. The upper and 
lower sidebands were self-calibrated separately and the so- 
lutions applied to the respective line data. 

The data reduction was performed using the MIRIAD 
package (Sault et al. 1995). A Briggs robust parameter of 
between and 2 was used for image construction and anal- 
ysis. For the purposes of presenting the data, the highest 



resolution with good signal-to-noise has been used (with 
a robust value of typically 0.25-0.5). Table 1 hsts the 
relevant source details, including the resolution of the im- 
ages presented. Flux densities in Jy may be converted to 
brightness temperatures in K assuming the following con- 
version factors for C^*0 and H2S respectively - 2.5 & 3.4 
(G9.62), 3.7 & 7.3 (G10.47), 3.9 & 4.3 (G29.96) and 6.4 & 
9.4 (G31.41). 

3. RESULTS 

In this section we describe the main observational re- 
sults, place them in context with previous observations 
and present the analysis alongside. C^^O and H2S spectra 
for each source are shown in Figure 1, clearly showing the 
high- velocity wings in the H2S spectra toward three of the 
four target regions. The wings arc most prominent toward 
G9.62 and G29.96, while the broadest lines are seen toward 
G10.47 and G31.41. Each source will now be described in 
more detail. 

3.1. G9.62+0.19 

The region around G9. 62-1-0. 19 is known to have nearly a 
dozen Hii regions (Cesaroni et al. 1994; Tcsti ct al. 2000), 
many of which are compact or ultracompact. Compo- 
nents D-G appear to lie within an elongated core, traced 
by CH3CN and continuum emission at 2.7 and 1.4 mm 
(Hofner ct al. 1996; Paper II). The brightest dust emission 
at 1.4 mm coincides with the centimeter source G9.62F, 
with a secondary maximum toward G9.62E. Testi et al. 
(2000) show that the centimeter fluxes of G9.62F arc con- 
sistent with an embedded B1-B1.5 star either exciting an 
ultracompact Hii region or driving an ionized stellar wind. 

The C^^O 2^1 emission towards G9.62 defines a single 
elongated core encompassing components F and D (Figure 
2a). Component E lies towards the edge of the integrated 
emission but does not stand out as a separate core. The 
emission shows a number of weaker extensions which do 
not appear to be associated with any of the known Hll 
regions. The overall distribution closely follows that of 
the C"'^*0 1^0 emission shown by Hofner et al. (1996). 
Channel maps (Figure 3a) reveal what may be the limb of 
a shell of gas defining the interface between the extended 
H II region B and the molecular cloud. A velocity gradient 
exists from north to south (blue to red) as seen previously 
by Hofner et al. (1996) in their 1^0 data. 

The C^^O spectra show self- absorption dips toward the 
center of the core. Comparing the total flux with the 
single-dish measurements of Hatchell et al. (1998b) show 
that our observations recover all of the C^^O emission to 
within the calibration uncertainties. Therefore, it is likely 
that these dips are genuine and not a result of missing flux 
in the map. The presence of the self- absorption makes it 
difficult to reliably estimate linewidths, but away from the 
core center the typical linewidth is ~4.0 km s~^ . There ap- 
pears to be a weak red wing on the C^^O spectrum toward 
component F, as also observed by Hofner et al. (1996). 
However, a map of the blue and redshifted emission does 
not show clear evidence for an outflow. 

The total integrated H2S emission in G9.62 also defines 
an elongated core (see Figure 2b), but it is clear that the 
denser gas traced by H2S is concentrated toward compo- 
nents E and F with two clear maxima seen toward these 
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sources. The emission also extends to cover components D 
and G but it is not possible to distinguish separate cores. 
The H2S emission peaks coincide with the 1.4-mm con- 
tinuum sources presented in Paper II. No H2S emission 
is detected toward components B, C, H or I. No emission 
was detected toward the extended near-infrared feature 
discovered by Persi et al. (2003). 

The H2S lines tend to be broader than the C^^O lines, 
indicating their formation in more turbulent regions closer 
to the center of the core, and arc broadest at the posi- 
tion of component F with a fuU-width-at-half-maximum 
(FWHM) of 6.8 kms-^ The FWHM linewidth at E is 3.1 
kms"""^, 5.2 kms"-'^ at G and 3.8 kms"-*^ toward component 
D. The spectrum toward F shows both blue and redshifted 
wing emission, and a map of the blue-shifted emission be- 
tween —6 and kms~^ reveals a well-defined lobe of emis- 
sion slightly displaced from the center of core F. The red- 
shifted gas (integrated between 8 and 13 kms"-'^) defines 
a counterpart to the blue lobe, and is centered on com- 
ponent F: see Figure 2c. Most of the H2O masers near 
G9.62F are redshifted, suggesting that they arise in the 
outflow. The outflow lobes do not show the same position 
angle nor do they peak at the same position as those seen 
in the HCO"*" data of Hofner et al. (2001) but they are in 
the same sense: blueshifted gas to the SW, redshifted to 
the NE. The H2S results therefore reinforce the interpre- 
tation that F is the hot core source and that it contains a 
young high-mass protostellar object. 

Like the C^*0, the H2S data also delineate a clear ve- 
locity gradient along the major axis of the core, also seen 
previously by Hofner et al. 1996), shown in a grayscale 
representation in Figure 2b (although not so obvious from 
the channel maps in Figure 3b). The H2S outflow emerges 
at an angle to the long axis of the core. The position angle 
of the outflow is 10.5 degrees east of north; that of the core 
is —13 degrees. (The position angle of the HCO+ outflow 
is 28 degrees.) The orientation of the blue and red lobes of 
the outflow deflnes a velocity gradient which is in the oppo- 
site sense to that seen along the major axis of the core, i.e. 
redshifted in the south and blueshifted to the north and 
west, further supporting our interpretation that the line 
wings arise in a bipolar outflow driven by an embedded 
source within F. Hofner et al. (2001) present arguments 
for the outflow being viewed pole-on. In this geometry it 
is unlikely that significant rotation can be detected. An 
alternative (and perhaps more likely) possibility is that 
the velocity gradient is due to the superposition of three 
cores (housing components D, F and G) which lie at dif- 
ferent velocities. Higher angular resolution would enable 
this hypothesis to be tested. 

While G9.62E is clearly detected in H2S emission, no 
line wings are observed which implies that either it has yet 
to form an outflow or it has evolved beyond the outflow 
stage. Given its well-defined H 11 region (stronger than F) 
and weaker dust emission than component F, it is likely 
that G9.62E is more evolved than G9.62F and has begun 
to clear away its parent core. 

3.2. GIO47+O.O3 

GlO.47-1-0.03 houses a group of four ultracompact Hll 
regions, three of which (A, Bl and B2) lie within a re- 
gion 2 arcsec across (C98). C98 argue that the majority 



of the NH3 (4,4) emission originates behind the three H 11 
regions, with component A being the least embedded. The 
ammonia peaks approximately 0.5 arcsec east of Bl. The 
fourth H II region, G10.47C, lies outside the area of ammo- 
nia emission. In this interpretation, the radio continuum 
sources are not deeply embedded within the hot core, in- 
stead lying close to the front surface. 

The C^^O emission (Figure 4a) defines a single com- 
pact core peaking on the Hll regions Bl and B2, with 
an extension to the west. No emission is dcitected toward 
component C. A number of features including the west- 
ward extension show up more clearly in the channel maps 
(Figure 5a). The westward extension covers the veloc- 
ity range 61.5 to 64.5 kms^^ (curving further north at 
64.5 kms~^), while between 64.5 and 69 kms^^ a corre- 
sponding feature exists in the south-east to east, with a 
well-defined tongue of emission extending to the east at 
69 kms~^. Such features could represent a bipolar out- 
flow, or a rotating (fragmented) core. We rule out outflow 
since the spectra to the east and west of center peak at 
different velocities, and thus the east and west extensions 
are not due to high- velocity line wings. Therefore it seems 
likely that these features define a gradient in the bulk of 
the core, and may represent smaller and less-massive indi- 
vidual cores. 

The C^^O spectra are single peaked across the entire 
region, and broadest at the position of the hot core. The 
lincwidths at the positions of the cast and west compo- 
nents are much narrower than toward the hot core (~3 
kms~^). Comparing our data with those of Hatchell et al. 
(1998b) again shows that we detect all of the C^*0 flux to 
within the calibration uncertainties. 

The H2S line is only detected in emission towards 
G10.47 and we do not see the same bow-shape seen in 
the ammonia of C98 caused by absorption against the H 11 
regions. The integrated H2S emission (Figure 4b) shows 
a single-peaked core elongated in a NW-SE direction at 
the lowest contours, but twists to a NE-SW direction in 
the highest contours, encompassing the UCH 11 regions, Bl 
and B2. This suggests that there are two separate cores 
containing Bl and B2 which we are unable to resolve. Like 
C98, we also do not detect any molecular line emission to- 
ward G10.47C. 

The brightest emission lies slightly south of Bl/2, and 
extends further south by ~2 arcsec, consistent with elon- 
gation in dust continuum (Paper II). The spectra toward 
each of the H 11 regions are flat- topped which could be due 
to either the presence of two velocity components which 
we are finable to resolve (perhaps from an expanding shell) 
or high optical depths in the H2S line. In our discussion 
below (Section 4.1) we conclude that the H2S is most likely 
to be optically thick. 

The H2S covers a velocity range from 58 to 74 kms~^ 
a similar range to the C^^O. Unlike the other sources in 
this sample, the H2S line does not show the same clear 
high-velocity wings (Figure 1). However, the H2S lines to- 
ward G10.47 are the broadest in our survey with a peak 
FWHM of 10.5 kms~^. From channel maps it is possible 
to distinguish two velocity components, although these do 
not correspond to any of the kinematic features of the 
C^^O data. The first component extends in the NE-SW 
direction and covers the range 61 kms~^ to 73 kms~^. A 
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second component extends NW-SE which traces material 
between 59 and 70 kms"'^. The NE-SW velocity compo- 
nent is also seen in NH3 (4,4) of C98. A secondary peak 
is observed at (3,-7) arcsec offset over the velocity range 
65 to 67 kms^^, coincident with one of the C^^O cores 
desrcibed above. 

While the H2S spectra do not show clear high- velocity 
line wings, the extreme red and blue shifted emission does 
show a slight north-south separation, although it is less 
than one beam diameter (Figure 4c). The center of this 
'outflow' lies ~ 1" north of the the 1.4-mm continuum 
peak and the UCHii regions Bl and B2. The fact that all 
the water masers are aligned on a north-south axis lends 
support to an outflow interpretation. The outflow inter- 
pretation gains further support from the fact that Olmi 
et al. (1996) detected a velocity gradient in ^^CO consis- 
tent with a north-south outflow (in the same sense as we 
observe in H2S), as well as an cast-west velocity gradient 
in CH3CN, i.e. perpendicular to the proposed outflow (al- 
though see § 5.2.1 below). The small offset between the 
lobes means that we have not resolved a flow direction, or 
that there is no outflow and we are be seeing the front and 
back surfaces of an expanding shell of gas centered on the 
H II regions. We therefore conclude that there is only weak 
evidence for either outflow or rotation toward G10.47. 

3.3. G29.96-0.02 

The ultracompact Hii region in G29.96— 0.02 shows a 
classic cometary shape (Wood & Churchwell 1989), which 
is also evident in our 2.7 and 1.4-mm continuum data (Pa- 
per II). C98 also observed G29.96 in ammonia (4,4) and 
detected a single peak coincident with a group of water 
masers, but offset from the H 11 region. Dust emission 
peaks at the position of the ammonia and water masers 
(Olmi et al. 2003; Paper II) confirming the location of the 
hot core. 

The C^^O emission toward G29.96 defines an extended 
core with multiple peaks (Figure 6a). There is general 
good agreement with the 1^0 data of Olmi et al. (2003), 
although they do not detect the westward ridge we see in 
our data. The C^^O emission shows a well-defined edge ex- 
tending from north-east to south-west, which is likely asso- 
ciated with the boundary between the molecular cloud and 
the ultracompact Hll region. This is seen more clearly in 
the channel maps. Figure 7a. This fact indicates that the 
hot core and the ultracompact Hii region are embedded 
within the same molecular cloud. The lone water maser 
~4.5" to the SW does not appear to be associated with 
any other peaks in either the C^^O or H2S emission. 

The C^^O emission peaks on the hot core, decreasing 
in intensity to the west, where at least two further cores 
can be discerned. These cores have not been resolved by 
any previous observations. There is a very slight velocity 
gradient along this ridge with the more distant core being 
marginally blue-shifted relative to the hot core (spectra 
peak at 96.8, 96.5 and 96.1 kms~^ from east to west). 

The C^^O lines are relatively narrow across much of the 
source, peaking at the hot core position. G29.96 has the 
smallest range of C^^O emission of our sample of sources, 
covering only 5 km s~^ (Figure 7a). The linewidths de- 
crease along the westward ridge away from the hot core 
(4.0, 3.5 and 3.2 kms'^ for the hot core, Wl and W2 re- 



spectively). Once again, there is good agreement between 
our C^^O flux and the single-dish measurement of Hatchell 
et al. (1998b). 

The integrated H2S emission shows a single maximum 
coincident with 1.4-mm contimmm peak (Figure 6b). 
Channel maps show that the H2S is extended over a region 
^ 15" in diameter and appears to wrap around the head of 
the UCHII region (Figure 7b), in a similar manner to the 
C'^^O. The H2S emission in G29.96 is the most extensive 
of all our four sources, and shows extensions to the west 
and south-west which may be further individual cores. 

The H2S spectra (Figure 1) show strong line wings ex- 
tending to ±10 kms""'^ across a wide area. Redshifted 
wings are seen in the south-east while blue-shifted wings 
extend to the north-west. Figure 6c shows the integrated 
red and blueshifted emission and clearly shows the spatial 
bipolarity of an outflow. The outflow center appears to lie 
closer to the ammonia peak of C98 than our 1.4-mm peak. 
This may indicate that there are multiple sources within 
the G29.96 hot core: higher resolution observations are 
clearly needed. The orientation of the blue lobe is along 
line of methanol masers seen by Minier, Booth & Conway 
(2002). The methanol masers are also blueshifted, increas- 
ingly blueshifted with distance from the outflow center. 

Figure 6b also shows an image of the first moment of 
the H2S emission (clipped at a level which excludes the 
wing emission), revealing a velocity gradient along a SW 
to NE axis, blue-shifted in the SW and redshifted toward 
the NE. There is evidence of this same gradient in the 
C^^O data as well. Figure 7a shows that the C^*0 peaks 
to the SW of the hot core at 95.5 kms"-'^ and to the NE at 
97.7 kms~^. It is tempting to interpret this as rotation of 
the hot core, especially since the gradient is observed to be 
perpendicular to the outfiow direction. However we regard 
this conclusion as somewhat tentative and clearly requires 
further observations of high-density tracers which arc not 
contaminated by outfiow. Under the assumption that we 
are detecting rotation the velocity gradient is ~11 kms"-^ 
pc~^, corresponding to a rotation period of 8.9 x 10^ yr. 

3.4. G3141+0.31 

The region associated with G31. 41-1-0. 31 houses an ul- 
tracompact Hii region along with a diffuse shell of ex- 
tended radio continuum emission (Wood & Churchwell 
1989). Hofncr & Churchwell (1996) detect a group of 8 
water masers, none of which coincide with the UGH 11 re- 
gion. However, the ammonia observations of C98 reveal a 
compact core of molecular gas clearly coincident with the 
water masers, but with no associated centimeter contin- 
uum emission. 

Our C^*0 emission shown in Figure 8a defines an elon- 
gated core, mostly oriented north-south but with a clear 
contribution from a diffuse, more extended component ori- 
ented NE-SW. Most notable is the fact that the C^^O does 
not peak at the hot core position; in fact the hot core is 
seen as a local minimum in the G^^O emission. A compari- 
son of our data and the JCMT measurement of Hatchell et 
al. (1998b) shows that we only detect of order 70% of the 
G^^O flux. It may be that the C^^O is not strongly peaked 
on the hot core position and is thus filtered out by our 
interferometric observations. Alternatively, the brightness 
temperature of the continuum emission may be sufficiently 
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high (at 27.5 K for the naturally-weighted data) that the 
C-'^^O is simply not detectable above this level. There are 
a number of peaks evident in the integrated map, although 
it is not certain if they represent independent cores given 
that some of the flux is missing from the image. 

Channel maps reveal the structure of the core. Figure 
8a shows the presence of a clear velocity gradient along the 
major axis of the C-'^^O core, blue-shifted in the south and 
red-shifted to the north. The individual peaks seen in the 
integrated intensity map are clearly seen in the channel 
maps. 

Figure 8b shows that the integrated H2S emission is cen- 
tered on the 1.4-mm and ammonia emission, with the wa- 
ter masers mostly lying along an axis parallel with the 
observed extension of the H2S core. The first moment im- 
age also shown in Figure 8b shows a weak velocity gradient 
along the core (in the same sense as the larger-scale C^^O 
gradient), which is tempting to interpret as rotation. The 
high-velocity wing emission (shown in Figure 8c) is pre- 
dominantly east-west in orientation. A position-velocity 
cut along the core shows the velocity gradient more clearly 
(sec § 5.2 below). We need further observations of com- 
plementary tracers to isolate one velocity component from 
the other. We note that the lower-resolution HCO"*" 1^0 
and SiO 2— >1 observations of Maxia et al. (2001) show the 
same distribution of red- and blueshifted high- velocity gas, 
supporting the outflow interpretation. The water masers 
appear to be associated with the outflow since many of the 
maser spots are red-shifted. 

Channel maps (shown in Figure 9b) show that the H2S 
emission is largely extended in the same NE-SW direction 
at all blue-shifted velocities. At the south end of the core, 
two secondary peaks show up at 91 to 93 kms~^, and it 
may be these components which bias the first moment im- 
age toward greater blue-shifted velocities in this regions. 
These may represent separate independent cores since our 
C^^O data show an extended maximum toward these po- 
sitions (and at the same velocity), as do the HCO"'" data 
of Maxia et al. (2001). The H2S spectra show varied line 
shapes across the core, with clear blue- and red-shifted 
wing emission toward some positions, and double-peaked 
line profiles toward the center of the core. However, as de- 
scribed above for G10.47, the double-peaked profiles may 
be due to self- absorption of an optically thick line. 

4. ANALYSIS: HOT CORE MASSES AND OUTFLOW 
PARAMETERS 

In this section we describe our analysis of the data in 
order to derive various physical parameters of the hot 
cores. We have made two-dimensional gaussian fits to 
the integrated-intensity images to derive source dimen- 
sions and intensities. In most cases a gaussian fit was 
a good representation of the observed emission. The de- 
convolved source dimensions were estimated assuming a 
gaussian geometry: these are the values listed in Table 
2. The mean linewidth was derived from gaussian fits to 
spectra toward each source (also given in Table 2). Next, 
virial masses were calculated assuming that each core can 
be approximated to a sphere of constant density. If the 
density falls off with radius then, for typical power-law in- 
dices, the virial masses will be smaller by up to 40 per cent 
(see e.g. MacLaren, Richardson & Wolfendale 1988). H2 



column and volume densities were calculated from these 
virial masses, again assuming constant density spheres. 
For C^^O we also calculated a total H2 mass based on the 
integrated intensity since the abundance is assumed to be 
well known (2xl0~^). These derived values are listed in 
Tables 3 and 4. 

4.1. Masses and column densities 

There is generally good agreement between the virial 
mass estimates from the H2S data and those derived from 
the C^^O. In the C^^O data, G29.96 has the lowest virial 
mass and G10.47 the largest. G10.47 also has the largest 
Il2S-derived virial mass, although now G9.62E has the low- 
est value. With the exception of G29.96, the C^^O virial 
masses are larger than the H2S masses, a result of the 
larger core size measured in the C^^O data. As expected, 
the H2S is tracing denser gas than the C^^O, and yields 
larger estimates of the H2 column density on account of 
the smaller source size. 

The molecular hydrogen column densities inferred from 
the H2S virial masses are exceedingly high. In particu- 
lar, G10.47 has the highest column density in our sample 
at 6 X lO^^cm"^, corresponding to a visual extinction of 
12 600 magnitudes. Even at mid-infrared wavelengths, the 
extinction remains significant. Using the values given by 
Mathis (1990), this H2 column density corresponds to an 
extinction of 270 300 magnitudes at 20 fim. The extinc- 
tion extrapolated to 1.4 mm is only 0.17 0.20 magnitudes, 
or essentially optically thin (t ^0.25). However, since the 
grain properties in the dense regions toward these hot cores 
can be different to those in the interstellar medium, the 
actual submillimetcr optical depth could be higher. The 
H2S virial masses imply a mean density in the four hot 
core sources of order 10''cm~^, in good agreement with 
values derived by Cesaroni et al. (1994). 

We have calculated the C^^O and H2S column density 
toward the hot core positions in each source (plus com- 
ponent E in G9.62), shown in Tables 3 and 4. We have 
assumed the emission to be optically thin, and that the 
level populations can be characterized by a single excita- 
tion temperature. The choice of excitation temperature 
is not constrained since we only have a single transition. 
For C^^O we have assumed values between 16 (equal to 
the energy of the J=2 level to obtain the minimum value: 
Macdonald et al. 1996) and 125 K (a reasonable upper 
limit based on ammonia and CH3CN estimates, e.g. Ce- 
saroni et al. 1994; Olmi et al. 1996). Using the H2 col- 
umn densities derived from the C^^O virial masses to es- 
timate an abundance of C^*0 yields values smaller than 
the canonical value of 2 x lO"'^ for G10.47 and G31.41, but 
in reasonable agreement for G9.62 and G29.96. Therefore, 
either the C^^O towards G10.47 and G31.41 is not opti- 
cally thin, or we are resolving out a significant portion 
of the C^^O emission. Since the agreement with previous 
single-dish observations is good for G10.47, we conclude 
that the C^^O is optically thick in this source. The C^^O 
and C^^O data of Hofner et al. (2000) show that the C^^O 
2^1 line has moderate optical depth (^1-2) in an 11" 
beam for all our sources. However, it is likely that the 
missing flux contributes to the poor agreement for G31.41. 
We have also calculated masses from the integrated inten- 
sity maps, assuming the standard C^^O abundance and 
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optically thin emission. Tabic 3 shows that there is good 
agreement between the masses calculated in this way and 
the virial masses for G9.62 and G29.96. The integrated 
masses are much lower than the virial masses for G10.47 
and G31.41, which is probably due to the effects discussed 
above. 

The observed H2S line brightnesses yield minimum ex- 
citation temperatures of 18 to 40 K (Figure 1). However, 
a strict lower limit for H2S can be obtained assuming that 
Tex = 2£'u/3fc = 56K (Macdonald ct al. 1996) where 
is the energy of the upper level, equivalent to 84 K in 
temperature units for the H2S 2(2,0)— >2(1,1) transition. 
However, the exact value is relatively unimportant as the 
derived H2S column density changes by only 50 per cent 
over the temperature range 30 to 250 K. Therefore, our 
choice of excitation temperature results in only a small 
uncertainty in our calculations, and lower than the uncer- 
tainty due to unknown optical depth. As a result, we only 
list a minimum value in Table 4. Examining the line in- 
tensities we see that for an excitation temperature of 56 
K, the lines in G10.47 and G31.41 are not optically thin 
(t ~ 1, and higher for lower excitation temperature). 

The H2S column density varies from 3 x 10"'^^ cm~^ to- 
ward G9.62E to 4.1 X 10^*^ cm'^ toward G10.47. If we as- 
sume that the virial masses calculated above are accurate 
representations of the mass of the region containing the 
H2S then we can estimate the abundance of H2S in each 
source. This yields column-avoragcd H2S abundances of 
1.2-6.8 X 10~^ relative to H2. Bearing in mind the fact 
that the virial masses are probably upper limits, and the 
fact that the H2S column density is a lower limit, these 
abundances should themselves be regarded as lower lim- 
its, although they are in good agreement with estimates 
from Hatchell et al. (1998b) and van der Tak et al. (2003). 

4.2. Comparison with previous estimates 

The first comparison we make is with our own BIMA 
continuum data presented in Paper H. To estimate the 
dust masses we assumed an dust emissivity of 0.018 
cm^g~^ at 1.4 mm (Ossenkopf & Henning 1994) and 
a temperature between 50 and 250 K. This yields dust 
masses in the range 17 to 94 Mq (G9.62F), 151 to 824 
M0 (G10.47), 36 to 198 Mq (G29.96) and 156 to 850 Mq 
(G31.41) with the larger values corresponding to the lower 
temperatures. In all cases, the C^**0 and H2S are consider- 
ably more extended than the dust emission, with FWHM 
dimensions typically twice as large as those of the dust 
cores. 

One the whole, there is reasonable agreement between 
the dust masses and the virial masses for both C-'^^O and 
H2S, at least for the lower dust temperatures. If the dust is 
indeed as hot as 250 K, however, then the discprcancy be- 
tween the dust and molecular line estimates can be large, 
ranging from a factor of 3 to more than an order of magni- 
tude difference. If the dust emission is optically thick, wc 
will underestimate the dust masses. Earlier we estimated 
a dust optical depth of 0.25 at 1.4 mm, which would lead 
to us underestimating the mass by only 10% or so. On 
the other hand, the H2S is clearly associated with outflow 
motions, which artifically raises the virial mass (due to 
increased linewidths), leading to the observed large dis- 
crepancy between the mass estimates. Using the (lower) 



mass estimates from the dust emission to derive H2 col- 
umn densities will lead us to derive higher values for the 
H2S abundance. Performing this calculation yields H2S 
abundances a factor of 3-5 higher than given in Table 4 
(again dependent on the assumed dust temperature). 

The H2S virial masses for both component G9.62E and 
G9.62F are in good agreement with the masses derived 
from Ci^O by Hofner et al. (1996), and our C^^O virial 
mass agrees well with their total mass for E and F. Since 
we measure a larger lincwidtli for oiu' H2S data (7.8 kms^^ 
compared with 4.4 kms~^ for ammonia), our estimate of 
the virial mass for F is correspondingly higher than the 
Ccsaroni et al. (1994) value although the sizes are in good 
agreement. Consequently, we also derive higher densities. 
Our larger linewidth is a reflection of the presence of the 
H2S in the outflow; ammonia is therefore not present or 
excited in detectable quantities in the outflow. 

Mass estimates for the other three sources vary consid- 
erably. Virial masses derived from the ammonia (4,4) line 
(Ccsaroni et al. 1994) tend to be ~5 times smaller than our 
values, largely due to the higher H2S linewidth and larger 
extent in C^^O. Tracers employed at 3 mm by Maxia et 
al. (2001) and Olmi ct al. (2003) yield estimates similar 
to what we derive. Maxia et al. (2001) use a smaller dust 
emissivity that we have employed; consequently we revise 
their estimates derived from 3-mm dust continuum down- 
wards by a factor of ~4 to give ~700 M© and 350 Mq 
for G29.96 and G31.41 respectively Olmi ct al. (2003) es- 
timate a mass of 1200 Mq for G29.96 from C^^O, larger 
than our estimates (Table 3, 4). 

Submillimeter dust continuum models constructed by 
HatcheU et al. (2000) for G10.47 and G31.41 have a con- 
stant density core at the center which has a mass of 3000 
Mq within a radius of 0.06 pc, somewhat larger than our 
estimates for a similar-sized region. 

The dimensions we derive for G9.62F are very similar 
to what Ccsaroni et al. (1994) determine from their am- 
monia data. In G10.47, G29.96 and G31.41 the size of the 
hot core derived from the H2S emission is larger than that 
derived from the ammonia observations of Cesaroni et al. 
(1994) and C98. C98 showed that the optically thicker 
main hyperfine lines of ammonia yielded a larger source 
size than the satellites. Since we find that the H2S emis- 
sion is probably optically thick, it seems plausible that the 
larger source size is due to a high optical depth in the H2S 
lines. Alternatively, our BIMA observations may be more 
sensitive to the distribution of H2S relative to the VLA 
observations of ammonia (Cesaroni et al. 1994). A third 
(and perhaps most likely) possibility is that the ammonia 
(4,4) line simply traces warmer gas than the H2S as the 
energy of the upper level is 84 K for the H2S while that of 
the NH3 (4,4) line is 201 K. 

4.3. Outflow parameters 

In Table 5 we list the derived properties of the outflows 
toward each source, including the mass of H2S-emitting 
gas. If we assume an abundance, then we can estimate 
the mass of swept-up material in the outflows from each 
source on the condition that the H2S is well- mixed with 
the molecular hydrogen. The H2S abundance is not well 
known in molecular clouds, and even less well known in 
outflows. In the previous section, we derived column- 
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averaged abundances of order 10^^, which are probably 
lower limits. In other work, Minh et al. (1991) derived 
values close to 10 ~* in a number of high- mass star- forming 
regions, and as high as 10~^ in Orion (Minh et al. 1990). 

For the gas in the outflow, we assume an abundance of 
10"'' relative to molecular hydrogen as it lies in the mid- 
dle of the range of values derived by Minh et al. (1990, 
1991), and 10~^ is used as an initial value for the gas- 
phase H2S abundances in the majority of chemical models 
(e.g. Rodgcrs & Charnley 2003, Wakelam et al. 2004). If 
the H2S abundance is as high as 10~^, then we derive mini- 
mum outflow masses of 0.4 Mq, values typical of relatively 
low-mass stars. On the other hand, if we assume smaller 
values then the mass estimates will increase to very high 
values (40 Mq per lobe) . Therefore, we regard our assump- 
tion of 10"'^ for the H2S abundance as providing reason- 
able lower-limits to the outflow parameters, although CO 
observations are desireable to determine these parameters 
with greater certainty. 

It appears that for all our sources several solar masses 
of material has been swept up, typically ~4 Mq in each 
lobe. Since we have neglected emission close to the line 
center, the true mass of each outflow is probably higher. 
Even so, given the mass estimates in Table 5, it is clear 
that these outflows are being driven by massive protostars 
since they compare favorably with a number of measure- 
ments of larger, more well-developed flows from massive 
YSOs (e.g. Gibb et al. 2003; Shepherd et al. 1998). 

If we assume that the outflow masses we derive are a 
sufliciently accurate reflection of the true masses, we can 
also estimate other properties of the outflows such as the 
momentum and energy. We assume that the characteristic 
outflow velocity is fout = Wmax/2 so that the momentum, 
p, is simply mVout, and the energy is Ek = Sim- 
ilarly, we assume the outflow 'size" (r-Qut) to be half the 
mean diameter of the outflow lobes. The dynamical 'age', 
tdyn; is calculated from rout/wout, with outflow force and 
mechanical luminosity calculated as p/tdyn and i^K/tdyn- 
Table 6 lists the values thus calculated. The average value 
for each outflow is given; the total momentum, energy, etc., 
will be a factor of two higher. All the values are consistent 
with the outflows being driven by relatively massive YSOs 
(c.f. Shepherd & Churchwell 1996), although the choice of 
H2S abundance clearly affects our flnal numbers. However, 
we wish to point out that since we are not sensitive to the 
highest-velocity gas (the velocity extent is quite low) the 
momentum and energy for each outflow will be seriously 
underestimated. High- J CO observations are desirable in 
order to better constrain the flow energetics. 

5. DISCUSSION 

H2S has only been previously observed in single-dish ob- 
servations of these sources. Since we are able to resolve the 
structure of the H2S emission we can now discuss our ob- 
servations in the context of physical and chemical models 
of hot cores. 

5.1. H2S in hot cores: distribution and abundance 

Hatchell et al. (1998a) estimate the H2S column density 
toward all four of our target sources from single-dish ob- 
servations of the same transition (and assuming the same 
excitation conditions) as we present here. These single- 



dish column densities are typically a factor of 20 lower 
than the column densities we derive, due to the bright- 
ness temperature being diluted by the 22-arcsec JCMT 
beam. Thus the filling factor for the JCMT observations 
is ^1/20, which predicts a source size of order 5", very 
similar to what we observe (Table 4). 

In the absence of any other input, chemical models of 
hot cores start with a gas-phase H2S abundance of around 
10"'' based on the observations of Minh et al. (1991), and 
more recently on the ISO results which place an upper 
limit of 10"'' on the solid H2S abundance (van Dishoeck 
& Blake 1996). For example, Rodgers & Charnley (2003) 
follow the evolution of chemical abundances in centrally- 
heated cores assuming two different models for the evo- 
lution of the physical parameters of the core. Rodgers & 
Charnley suggest that their static model is more appropri- 
ate for massive cores, and show that the abundance of H2S 
is initially high and decreases rapidly with time. At early 
times (< 10^ yr) H2S appears to have a uniform abun- 
dance across the core, falling off at an outer radius which 
increases with time, reaching 0.01 pc at lO"* yr (the outer 
radius of their model is 0.03 pc). The typical abundance 
for II2S in these models is a few x 10^*^ to 10^^, falling off 
rapidly with radius at late times (although it rises again 
in the outer regions of the core). Inspection of their pre- 
dicted column densities shows that only at late times (at 
10^ yr and beyond) does the H2S column density fall be- 
low 10^^ cm^^. Simultaneously, the size of the hot core 
in ammonia is at least a factor of two smaller, purely as a 
result of the chemistry. This appears to be another factor 
which can account for the fact that these hot cores appear 
to be smaller in the ammonia images of C98 than in our 
H2S. 

We note that the hot core models predict H2S abun- 
dances at least an order of magnitude greater than the 
values we calculate (Table 4). As discussed elsewhere in 
this paper, it is likely that a combination of high optical 
depth and an overestimated mass from the virial theo- 
rem can account for this discrepancy. Since ammonia is 
well studied on similar scales in hot cores, we can exam- 
ine the NH3/H2S abundance ratio. Using the ammonia 
column densities given by Cesaroni et al. (1994) we derive 
ratios [H2S/NH3]=2-8x 10"^. However, in the models of 
Rodgers & Charnley (2003), this ratio never falls below 
unity except at late times (10^ yr). The simplest expla- 
nation for the observed ratios is that the II2S lines have 
significant optical depth leading to vastly underestimated 
H2S column densities. 

5.2. Hot core dynamics: outflow and rotation? 

The data we present here appear to trace two kinematic 
components. The first is outflow traced by the H2S. The 
second is a velocity gradient more or less perpendicular to 
the outflow and is probably rotation, seen in both C^^O 
and II2S. The outflow evidence is clear: spatially bipolar 
lobes of red- and blueshifted gas. The rotation interpre- 
tation is ambiguous, especially in light of the likely frag- 
mented nature of the G9.62 core and the more obvious out- 
flow signature in the II2S. However, given the good agree- 
ment between the C^^O and H2S gradients across each 
core, the rotation interpretation seems more tenable. In 
Figures 10 to 13 we show selected position-velocity (PV) 
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cuts across each of the H2S data cubes for each of the 
hot cores (along the fong axis of the core and along the 
outflow axis where defined). In three of the four sources 
clear velocity gradients are present. Only G10.47 (Figure 
11) does not show a velocity gradient. The observed ve- 
locity gradients tend to be linear along the long axis of 
the core, although they typically only extend for a cou- 
ple of beamwidths. (However, we note that it is possible 
to distinguish velocity gradients on scales smaller than a 
beam: e.g. Olmi ct al. 1996.) It is clear that the veloc- 
ity gradients do not directly show evidence for Keplerian 
rotation; indeed given the massive cores within which the 
central stars are embedded such a result would be sur- 
prising since a Keplerian rotation law only strictly applies 
for a dominant point mass (although it may be true that 
the velocity field is locally Keplerian). Furthermore, if the 
observed emission is optically thick (as concluded above) 
then regardless of the intrinsic velocity field, the observed 
velocity gradient will be linear. 

However, even if the emission is not optically thick, our 
finite beamsize prevents us from directly observing a Ke- 
plerian rotation curve, instead smoothing out the central 
region into a linear gradient (sec Richer & Padman 1991). 
Therefore, we caution against over-interpreting these PV 
data, especially since the emission closest to the center is 
heavily contaminated by the high-velocity emission from 
the outflow. At best we can only conclude that in the cases 
of G9.62, G29.96 and G31.41 our C^^O and H2S data show 
evidence for velocity gradients along the long axis of the 
hot core, but not what form the velocity law takes. 

If we assume the linear velocity gradicmts arc due to uni- 
form rotation then we can estimate rotation periods, and 
then the masses required to bind the motions. For each 
of the three sources with clear velocity gradients (G9.62F, 
G29.96 and G31.41), the rotation periods for each core is 
similar at 1 2x10"'' yr. Using the core dimensions shown 
in Table 4 these periods require binding masses of only 15, 
60 and 154 (for G9.62F, G29.96 and G31.41 respec- 
tively). These masses are all less than both the virial and 
dust-derived masses and so that the cores are rotation- 
ally bound, and may even be still collapsing. The veloc- 
ity gradients along the outflow axes are larger and imply 
binding masses of 810, 2800 and 2360 Mq respectively, all 
much larger than the virial masses given in Tables 3 and 
4. Therefore the motions along the outflow axes are not 
bound to the hot cores. 

An interesting feature of the G31.41 PV diagram is that 
the H2S emission farthest from the center does not lie at 
the rest velocity. This indicates that the H2S is not very 
abundant at these distances from the hot core and falls 
below our detection threshold. From Figure 13 we can 
therefore determine an outer radius of the hot core of ~5 
arcsec or 0.19 pc (5.9 x 10^^ cm). 

5.2.1. An outflow origin for methyl cyanide? 

Comparing our results with those of Beltran et al. (2004) 
for G31.41, we find that our outflow direction is closely 
aligned with the long axis of their methyl cyanide emis- 
sion, which they interpret as a rotating disk or toroid. 
Their interpretation was a natural conclusion given the 
example of IRAS 20126-1-4104 (Cesaroni et al. 1999) which 
shows a clear disk-like component for the CH3CN emission. 



This ambiguity highlights the kind of problems faced when 
studying hot cores. However, we believe that our C-'^^O 
and H2S data display clear evidence for both outflow and 
(in two cases) velocity gradients perpendicular (or close 
to) to the axis of the outflow. The simplest interpreta- 
tion for these gradients is rotation of the core. The distri- 
bution of our H2S (and 1.4-mm continuum) data closely 
matches that of the NH3 (4,4) of C98. This leads to the 
interesting conclusion that Beltran et al. (2004) have actu- 
ally detected CH3CN in the outflow from G31.41. Further 
support comes from the CH3CN observations of G29.96 by 
Olmi et al. (2003) which are also consistent with emission 
from an outflow since the velocity gradient is extremely 
high (higher than we detect with our H2S data, thus re- 
quiring an even higher binding mass) and oriented mostly 
along the same direction as our outflow. It is interesting 
to compare the methyl cyanide and H2S results. H2S is 
clearly not confined to the outflow. On the other hand, 
it seems as if the CH3CN in these sources is dominated 
by an outflow component, suggesting that it is perhaps 
shock chemistry that is playing a greater role in the for- 
mation of methyl cyanide. Clearly further high-resolution 
observations are necessary in order to explore this further. 

5.3. Embedded sources in hot cores 

The primary result from these observations is the con- 
firmation that outflows emanate from some hot cores, in- 
cluding all four of the targets in our sample. Further- 
more the images of the high-velocity emission show that 
these outflows appear to be coUimated; in other words, un- 
like many ultracompact H II regions, they are not spherical 
shells driven by the expansion of a symmetric stellar wind. 
This has important implications for the formation of mas- 
sive stars and implies the presence of a collimation mech- 
anism. In the case of low-mass stars, this is now widely 
believed to be a rotating circumstellar disk (Ouycd, Clarke 
& Pudritz 2003; Konigi & Pudritz 2000). Explanations of 
outflows from more luminous protostars invoke a wider- 
angle contribution from a stellar wind (Shepherd et al. 
1998; Richer et al. 2000), although a disk is still believed 
to play an important role in at least initially collimating 
the outflow. 

Clearly we cannot determine whether circumstellar disks 
exist around these outflow sources with the current ob- 
servations. However, we can examine whether the ambi- 
ent medium has sufficient turbulent pressure to act as a 
collimating agent for the flows. Assuming a single driv- 
ing source we may equate the momentum flux from the 
wind/outflow (assumed to be initially isotropic) with the 
turbulent back-pressure of the ambient medium: 



where p is the ambient density, afj-, is the 1-dimensional ve- 
locity dispersion {— -\/8 In 2 Aw, where Av is the FWHM), 
p is the outflow force and r is the radius at which the forces 
are equal (the collimation radius). Using values tabulated 
in Tables 4 and 5 we derive collimation radii of 350 to 950 
AU. At the distances to our sources, these radii correspond 
to angular scales of 0.06"to 0.13", and are thus unresolv- 
able with current facilities. If multiple sources contribute 
to the flow then the expression above can be applied for 
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each source assuming a momentum flux equal to the ap- 
propriate fraction of the total. This will lead to smaller 
coUimation radii and our conclusion is unchanged. There- 
fore, it seems likely that the ambient medium is playing a 
significant role in producing the bipolar appearance of the 
outflows from these hot cores. The fact that outflows are 
observed clearly points to a non-spherical geometry. 

Another implication from our observations is that the 
outflow must begin at a stage prior to the formation of an 
ultracompact Hll region since the hot cores in G29.96 and 
G31.41 have no radio sources within them. It could be 
argued that G10.47 provides evidence counter to this as- 
sertion, but we note that C98 conclude the H ii regions in 
G10.47 are located toward the front of the hot core, and are 
therefore not deeply embedded within the hot core itself. 
However, G10.47 is our weakest example of an outflow, 
and the hot core may only be a hot, dense remnant of the 
formation of the stars exciting the Hii regions, as Watt & 
Mmidy (1999) concluded for G34.26. Unfortunately it is 
not possible to assess the likelihood of this interpretation 
with our data and further observations are desirable. 

G9.62F has a radio source within it, but it is very weak, 
and probably represents an ionized stellar wind rather 
than an Hii region (Testi et al. 2000). The Orion hot 
core is probably excited by the radio source I (Menten & 
Reid 1995), which has a 22-GHz fltix density of 4.8 mJy 
(Plambeck et al. 1995). If source I were located at a dis- 
tance typical of the sources in this paper (^6.5 kpc) then it 
would be (currently) undetectable at 22 GHz (the 22-GHz 
flux density would only be 30 /LiJy). Equally, BN would 
also not be detectable. Therefore it is possible that equiv- 
alent sources are located within the hot cores studied here. 
It should be noted that the centimeter radio emission in 
both BN and source I is from an ionized stellar wind rather 
than an ultracompact H II region. 

Could it be that there arc Hii regions in G29.96 and 
G31.41 but they are not detectable due to their greater 
distance? We can test this by scaling the fluxes of the H ii 
regions in G10.47, and using the values given in C98 for 
the 1.3-cm flux, we flnd that such Hii regions would be 
detectable at the distance of G29.96 and G31.41. There- 
fore the non-detection of Hii regions toward these hot 
cores suggests that ultracompact H ii regions have not yet 
formed within them, although stellar wind sources may 
be present. Since the outflow phenomenon is intimately 
linked to accretion, our observations are also compatible 
with the hypothesis that Hii regions can be quenched by 
a strong accretion flow (Walmsley 1995). 

Another interesting point to note is that in each case, 
there appears to be only a single bipolar outflow detected. 
(If the observed bipolar outflow is the superposition of 
multiple flows, then they clearly share a common orienta- 
tion.) This is surprising given the clustered nature of mas- 
sive star formation: it seems more likely that we woiild not 
detect an orderly bipolar flow from such a crowded envi- 
ronment. Indeed, even in binary systems the outflow orien- 
tations can be distinctly different: the low-mass systems in 
HHlll and NGC 1333-IRAS2 both show pairs of outflows 
with very different position angles (Cernicharo & Reipurth 
1996; Knee & Sandell 2000), while IRAS 16293-2422 ex- 
hibits a quadrapolar flow (Walker, Carlstrom & Bieging 
1993). The observation of a single flow has important im- 



plications for the formation of massive stars. One pos- 
siblity is that the most massive and energetic flow (pre- 
sumably from the most massive YSO) dominates the en- 
ergetics, in a similar fashion to how infrared and radio 
continuum observations are dominated by the most lumi- 
nous source. Alternatively, the outflow stage could be so 
short that at any one time, we only catch one source in 
the act. An extension of this interpretation is that we have 
detected these outflows in H2S, a molecule which is itself 
comparatively short-lived in hot cores (Charnlcy 1997). It 
is possible that the short outflow lifetime coupled with the 
short chemical lifetime for H2S reinforce each other to pro- 
duce a single outflow detection. Of course it is highly likely 
that multiple sources are present within these hot cores, 
each potentially with its own outflow. High-resolution ob- 
servations of CO, whose abudance is well-known and sta- 
ble, are needed to search for multiple outflows. 

As always it is instructive to compare these results with 
the nearest hot core in Orion. The outflow from the Orion 
hot core is well-defined in CO (although the coUimation 
is not very high; e.g. Chernin & Wright 1995) but has 
not been probed in H2S. The CO 4^3 data of Schulz et 
al. (1995) show a bipolar outflow with outflow lobes sep- 
arated by ~20". If such a flow were moved ten to twenty 
times further away then it is the CO flow would only be 
barely resolvable with current interferometers. 

5.4. Scattered light at mid-IR wavelengths? 

An interesting feature of our G29.96 results is that the 
18 /xm peak of De Buizer et al. (2002) is offset to the NW 
of the 1.4- mm continuum peak (Figure 6c). It lies within 
the blue lobe close to the position of the methanol masers 
observed by Minier et al. (2002). This leads to the hypoth- 
esis that even at 18 /im we are seeing scattered light. If 
this is the case, then we expect the 10 /im emission to peak 
further from the outflow center. This does not appear to 
be the case but it is difficult to be sure since the hot core 
is weaker relative to the UGH 11 region at 10 /zm. This in- 
terpretation is tentative at best since the warm dust in the 
H II region dominates the emission at these mid-IR wave- 
lengths and makes the isolation of the hot-core emission 
non-trivial. However, it is useful to examine our data for 
further support. 

In Section 4.2 above, we estimated that the extinction 
toward the hot cores is significant even in the mid-infrared 
(450 to 520 magnitudes at 20 /im). This, and the extra 
extinction caused by the 18-/im silicate absorption feature 
(e.g. Mathis 1990), would result in the hot core emission 
to be strongly self-absorbed at 18-/xm. Therefore the only 
way that light can escape is through a cavity. An outflow 
provides a natural explanation for the presence of a cav- 
ity, and so we conclude that the evidence strongly suggests 
that the 18-/im source detected by De Buizer et al. (2002) 
is in fact scattered light. It is interesting to note the re- 
cent conclusion of McCabe, Duchene & Chez (2003) that 
the 12-/zm emission toward the T Tauri star HK Tau B 
appears to be scattered light. 

Hofner et al. (2001) argue that the outflow in G9.62F is 
observed almost pole-on since weak near-infrared emission 
is detected. If our scattered-light interpretation is correct, 
then G9.62F should be easily detected at 18-/im and ob- 
served coincident with the hot core, since the dust column 
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to the central source will be much less. Similar observa- 
tions of G10.47 and G31.41 arc desirable to test whether 
these sources are visible at 18 /um. 

6. CONCLUSIONS: COLLIMATED OUTFLOWS FROM 
ROTATING HOT CORES 

We have presented high-resolution observations made 
with the BIMA millimetre array of the 2(2,0)^2(1,1) tran- 
sition of H2S and 2^1 transition of C^^O toward a sample 
of four hot cores: G9. 62+0.19, G10.47-h0.03, G29.96-0.02 
and G31.41-I-0.31. In three of the four cases (G9.62F, 
G29.96 and G31.41), high-velocity H2S emission is ob- 
served which is distributed in a spatially-bipolar fashion, 
in a direction offset from the major axis of the hot core (as 
defined by the integrated intensity). In addition, linear ve- 
locity gradients are observed in both C^*0 and II2S along 
the major axis of these three cores. Only in GIO. 47-1-0. 03 
are the data ambiguous, showing neither clear rotation nor 
outflow signatures. 

Although the H2S abundance is not well known, we have 
made reasonable assumptions to derive the parameters of 
these outflows. These estimates show that the outflows are 
indeed driven by massive protostars despite being subject 



to large uncertainties in the II2S abundance. 

The lack of centimeter- wave radio emission from within 
two of the four hot cores shows that outflows begin at a 
phase prior to the generation of an ultracompact H 11 re- 
gion. Furthermore, our data show only one outflow toward 
each hot core. This implies that the observed outflow in 
H2S is dominated by a single source, either by the most 
massive or by the source which we just happen to catch in 
this particular phase of evolution. 

We interpret these data as indicating the presence of a 
massive protostar embedded within a hot core, which may 
be rotating, and driving a molecular outflow. The bipolar 
nature implies that these flows are coUimated, and we have 
shown that the ambient material is capable of channeling 
flows with the properties derived here and, if this is the pri- 
mary coUimation mechanism, then the mass distribution 
is not spherically symmetric. 
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Table 1 



Source details, coordinates for pointing centers and beam parameters. 



Source 


a(J2000)'' 


6(J2000)'' 


d 


fLSR 


Lhoi 


H2S 




















Beam 


PA 


Beam 


PA 








(kpc) 


(kms^'-) 


(Lq) 


(arcsec^ ) 


(deg) 


(arcsec^) 


(deg) 


G9. 62+0.19 


18:06:14.812 


-20:31:39.390 


5.7 


4.4 


4.4(5) 


3.4x1.7 


2.5 


4.5x2.3 


6.7 


GlO.47+0.03 


18:08:38.280 


-19:51:r50.000 


5.8 


67.8 


5.0(5) 


2.8x1.3 


9.2 


3.5x1.9 


8.6 


G29.96-0.02 


18:46:03.775 


-02:39:21.884 


7.4 


97.6 


1.4(6) 


2.0x1.3 


7.0 


3.5x1.9 


8.4 


G31.41+0.31 


18:47:34.320 


-01:12:45.800 


7.9 


97.4 


1.7(5) 


1.8x1.2 


24.8 


2.4x1.7 


24.5 



Note. — Bolometric luminosities are taJsen from Cesaroni et al. (1994) and Hatchell et al. (2000). The notation a(b) 
represents a x lO''. 

''In hours, minutes and seconds. 

''In degrees, minutes of arc and seconds of arc. 



Table 2 

Observed properties of each source C^^O and H2S. 



Source 


FWHM Dimensions'* 
(arcsec^) (pc^) 


PA 

(deg) 


(kms-i) 


H2S 

FWHM Dimensions^ 
(arcsec^) (pc^) 


PA 

(deg) 


(kms-i) 


G9. 62+0.19 


11.4x7.0 


0.32x0.19 


-33 


4.2±1.2 










G9.62+0.19E 










2.9x2.0 


0.08x0.06 


-29 


4.1+0.9 


G9.62+0.19F 










4.6x2.2 


0.12x0.06 


-24 


5.3+1.3 


GlO.47+0.03 


6.4x5.2 


0.18x0.15 


-42 


7.4±1.8 


5.5x3.6 


0.15x0.10 


-31 


8.2+1.7 


G29.96-0.02 


9.1x6.8 


0.33x0.24 


-70 


3.2±0.5 


4.3x3.2 


0.16x0.12 


-18 


5.2+1.9 


G31.41+0.31 


12.7x6.1 


0.49x0.23 


+22 


4.5±1.6 


3.5x2.5 


0.14x0.10 


+26 


6.6+1.8 



Note. — The notation a{b) represents a x lO''. 

^Deconvolved dimensions obtained from a two-dimensional gaussian fit to the integrated intensity 
image. The uncertainty is typically less than 0.25", or less than 0.01 pc at the distance to all sources. 

''Mean linewidth within the half-power contour. Uncertainty represents dispersion in values. 
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Table 3 



Parameters derived from C O data. 



Source 


(Mq) 


(cm ^) 


(cm~^) 


(cm~^) 


(Mq) 


G9.62+0.19 
GlO.47+0.03 

G29.96-0.02 
G31.41+0.31 


463±137 
920±252 
301± 52 
723±261 


7.7±2.3(23) 
3.8±1.0(24) 

4.0±0.7(23) 
6.5±2.3(23) 


1.0±0.3(6) 
7.7±2.1(6) 

4.6±0.8(5) 
6.2±2.2(6) 


1.2±0.6(17) 
1.8±0.8(17) 
9.7±4.5(16) 
1.5±0.7(16) 


303±141 
172± 80 
346±161 
264±123 



Note. — The notation a{b) represents a x 10^. 

''Calculated assuming a uniform sphere. 

Calculated assuming optically thin emission and an excitation temperature of 
75 K. Uncertainty reflects excitation temperatures in the range 16 to 125 K. 

'^At the hot core position 

■^Mass calculated within the half-power contour of the integrated C^^O emis- 
sion, assuming a C^^O abundance of 2x10"''. 



Table 4 

Parameters derived from H2S data. 



Source 












s 




(Mo) 


(cm-2) 


(cm-3) 


(cm-2) 






G9.62+0.19E 


124± 45 


2.6±1.0(24) 


1.2±0.4(7) 


>3.0(15) 


>1.2(- 


-9) 


G9.62+0.19F 


248± 85 


3.7±1.3(24) 


1.4±0.5(7) 


>1.2(16) 


>3.3(- 


-9) 


G10.47-h0.03 


861±228 


6.0±1.6(24) 


1.6±0.4(7) 


>4.1(16) 


>6.8(- 


-9) 


G29.96-0.02 


392±154 


2.2±0.8(24) 


5.0±2.0(6) 


>1.3(16) 


>6.1(- 


-9) 


G31.41-F0.31 


540±173 


4.1±1.3(24) 


1.1±0.4(7) 


>2.7(16) 


>6.7(- 


-9) 



Note. — The notation a{b) represents a x 10^. 
''Calculated assuming a uniform sphere. 

''Calculated at the peak position assuming optically thin emission and an 
excitation temperature of 56 K. 
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Table 5 
Outflow properties. 



Source 




Red lobe 






Blue lobe 






Sizc^ 


Velocity^' 




Size^ 


Velocity^ 






(pc2) 


(kms-i) 


(Mo) 


(pc2) 


(kms~^) 


(Mo) 


G9.62+0.19F 


0.12x0.07 


+6.5 


4.4(-7) 


0.14x0.08 


-8.5 


4.0(-7) 


GlO.47+0.03 


0.13x0.07 


+7.5 


3.5(-7) 


0.14x0.08 


-10.0 


5.4(-7) 


G29. 96-0.02 


0.13x0.08 


+7.5 


4.5(-7) 


0.15x0.09 


-12.2 


4.0(-7) 


G31.41+0.31 


0.11x0.07 


+11.9 


5.0(-7) 


0.07x0.04 


-13.8 


3.7(-7) 



Note. — The notation a{b) represents a x 10^. 

'^Deconvolved FWHM dimensions derived from a two-dimensional gaussian fit. 
''The maximum flow velocity relative to the source systemic velocity. 



Table 6 

Estimated outflow ENEROETicsf 



Source 




■Uout 


'dyu 


M 


M 


P 




P 




Lmoch 




(pc) 


(kms-^) 


(yr) 


(Mo) 


(Moyr-i) 


(Mq kms"^) 


(erg) 


(Mq kms" 


(L©) 


G9.62+0.19F 


0.10 


3.8 


2.6(4) 


4.2 


1.6(-4) 


16 


6.1(44) 


6.2(- 


-4) 


0.2 


GlO.47+0.03 


0.10 


4.4 


2.2(4) 


4.4 


2.0(-4) 


19 


8.5(44) 


8.8(- 


-4) 


0.3 


G29.96-0.02 


0.11 


4.9 


2.2(4) 


4.2 


1.9(-4) 


21 


1.0(45) 


9.4(- 


-4) 


0.4 


G31.41+0.31 


0.07 


6.4 


1.1(4) 


4.4 


4.0(-4) 


28 


1.8(45) 


2.6(- 


-3) 


1.3 



Note. — The notation a(b) represents a X lO''. 

Note. — Values listed are the average of the those for each lobe of the outflow. Totals will be a factor of two greater. 
^An H2S abundance of 10~^ relative to H2 has been assumed. 
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Velocity (km/s) Velocity (km/s) Velocity (km/s) Velocity (km/s) Velocity (km/s) 

Fig. 1. — C^^O and H2S Spectra toward the peak positions in each source. The vertical scale for each transition is the same in each panel 
to show the relative strength of the line toward e3x:h source. The dotted C'^^O profile for G31.41 represents the spectrum from the northern 
peaJi. 
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a) 



G9.62+0.19 
C^^O 2-1 integrated intensity 



□ HIT regions 

Water masers 
A Methanol maser 
• 1.4 mm peaks 




^\ G9.62+0.19 

HgS moment (contours) & 1 (gray) 
□ HII regions 
+ Water masers 
A Methanol maser 
• 1.4 mm peaks 




■4f- IR feature 



□ 



I 



5 0-5 
RA offaet (arcsec) 



^■j G9.62+0.19 

HgS blue and redshiftsd line winga 

□ HII regions 

+ Water masers -f^ 

A Methanol maser ^ LJ 

• 1.4 mm peaks 




□ 



RA offset (arcsec) 

Fic;. 2. — a) C^^O 2—^1 integrated intensity towards G9.62. Hll regions are marked with open squares and water masers from 
Hofner & Churchwell (1996) are marked by plus signs. The IR feature of Persi et al. (2003) is marked by a plus sign with a small 
filled circle. The size of the open squares reflects the centimeter-wave flux density: a larger square represents a greater fiux density, 
b) Contours of the zcroth moment (clipped at 2.0 Jybeam"^) superimposed on a color-scale representation of the first moment 
(clipped at 2.5 Jybeam^^). Contours begin at 6 Jy beam~^ kms"'^ and increase in steps of 4 Jybeam"^ kms~^. The grayscale 
scale varies linearly between 1.4 kms~^ (black) and 4.2 kms~^ (white). The beam is marked by the hatched ellipse in the bottom 
left-hand corner of the image, c) Superposition of blue-shifted (solid contours) and red-shifted H2S emission (dashed contours). 
Average emission over a 4 kms^^ -wide range centered at —4 kms~^ (blue) and 10 kms~^ (red). Contours increase in steps of 
0.36 Jybeam"^ starting at 0.54 Jybeam"^. 
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Fig. 3. — G9.62 channel maps of a) C^*0 and b) H2S emission from —5 to 11.5 kms~^ in steps of 1.5 kms~^. Contours begin at 1.2/0.45 
Jy beam^^ and increase in steps of 0.8/0.30 Jy beam^^ for C^*0 and H2S respectively. Hll regions are marked by open squares. The beam 
is marked by the ellipse in the bottom left-hand corner of the image. 
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a) GlO.47+0.03 

C^^O 3-1 integrated intensity 

n HII regions 
+ Water masers 




HA offset (a: 



b) 



GlO.47+0.03 
moment (contours) Sc 1 (gray) 

C 



□ HII regions 
-|-Water masers 
• 1.4 mm peak 




I 



5 0-5 
RA offset (arcsec) 



GlO.47+0.03 
HgS blue and redshifted line wings 




RA offset (arcsec) 

Fig. 4. — a) G10.47 moment zero contours (clipped at a flux density of 1.0 Jybeam"^) superimposed on a grayscale of the first moment 
(clipped at 2.0 Jybeam^-*^ to isolate the line core). Contours begin at 10 Jybcam^^ kms^-' increasing in steps of 10 Jybeam^-' kms^-*^. The 
peak is 78 Jybeam^-*^ kms^^. Ultracompact Hll regions are marked by open squares and water masers by plus signs. The beam is marked by 
the ellipse in the bottom left-hand corner of the image, b) G10.47 high velocity emission averaged over the velocity ranges 74 to 76 kms^^ 
(red) and 58 to 60 kms~^ (blue). Redshifted gas is shown by the dashed contours and blueshifted gas by the solid contours. Contour levels 
are the same for each lobe, beginning at 0.54 Jy beam~^ in increasing in steps of 0.36 Jybeam"^. 
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Fig. 5. — G10.47 channel maps of a) C^*0 and b) H2S emission between 57 and 73.5 kms ^. Contours begin at 0.6/0.75 Jybeam 
increasing in steps of 0.4/0.5 Jybeam^^ for C^®0 and H2S respectively. Symbols are as in Figure 4. The beam is marked by the ellipse 
the bottom left-hand corner of the image. 
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a) 



G29.96-0.0a 
C^^O 2-1 integrated intensity 




Water masers 
A Methanol maser 
• 1.4 mm peak 
18 fiTci peak 



b) 



5 0-5 
HA offset (arcsec) 

— 1 1 1 1 1 1 1 1 1 

Gag.g6-o.02 

moment (contours) He 1 (gray) 



□ Hn regions 
-|- Water masers 
A Methanol maser 
• 1.4 mm peak 
>|£ 18 /mi peak 




RA offset (arcsec) 



c) 



G28.86-0.02 
HgS blue and redshifted line winga 




RA offset (arcsec) 

Fic;. 6. a) Contours of zcroth moment superimposed on the first moment grayscale image (clipped at a flux density of 2 Jy bcani^-' to 
isolate the line core). Note the SW-NE gradient in the first moment data perpendicular to the direction of the outfiow lobes, suggestive 
of a rotating core, b) Red-shifted (dashed contours) and blue-shifted (solid contours) H2S emission averaged over a 4-kms^^-widc interval 
centered at 90 kms^^ (blue) and 106 kms^^ (red). The 18-/im peak is almost coincident with the methanol maser (De Buizer et al. 2002) 
which lies within the blue lobe of the outflow. Contours for each lobe begin at 0.1 Jybeam"^ and increase in steps of 0.1 Jybeam"^. The 
beam is marked by the ellipse in the bottom left-hand corner of the image. 
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64.0 b) G29.96 - HgS 

□ 4 

+ 

;Q 


86.0 

+ 

;Q 


88.0 

+ 

;Q 


90.0 

o 

+ 

;Q., 


93.0 

Pi 1 k 1 J 1 1 1 1 1 1 1 -1 1 1 1 1 1 4 


94.0 

;Q,..,^ 


96.0 


98. g 

,Q 


100.0 ' 

jQ,,;' 


103.0 

D 


104.0 

.0 


106.0 

1 

- 

+ 

, , , 



5 0-5 -10 5 

RA offset (arcsec; JSOOO) SA offset (i 



-5 -10 5 -5 -10 

JSOOO) RA offset (arcsec; JSOOO) 



5 -5 -10 5 -5 -10 5 -5 -10 

RA offset (arcsec; JSOOO) RA offset (arcsec; J8000) RA offset (arcsec; JSOOO) 



Fig. 7.— G29.96 channel maps of H2S emission from 90 to 102.1 kms-^ (C^^O) and 84 to 106 kms-^ (H2S) in steps of 1.1/2 
kms~^(C^^O/H2S resepctively). Contours begin at 1.2/0.3 Jybeam^^ and increase in steps of 0.8/0.2 Jybeam"^ for C^^O and H2S 
respecetively. The 1.4^mm continuum peak is marked by a filled circle, the position of the UCHii region is marked by an open square, water 
masers axe marked by plus signs, the methanol maser group of Minier et al. (2002) are marked by an open triangle. The beam is marked by 
the ellipse in the bottom left-hand corner of the image. 
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Fic;. 8. — a) Moment zero (integrated intensity, clipped at 0.5 Jybeam^^) contours superimposed on grayscale of the first moment 
image (clipped at 1.5 Jybeam^-"^). Contours start at 4 Jy beam~^ kms^-"^, increasing in steps of 5 Jybeam^^ kms^^. The peak is 41 
Jy beam^^ kms^^ . Grayscale is linear from 92.5 (blue) to 98.9 kms^^ (red) as shown on the sidebar, b) G31.41 redshifted contours (red) 
and blueshifted contours (blue) averaged over a 4-kms~^-wide interval centered at 89 kms~^ (blue) and 105 kms^^ (red). For each lobe, 
the contours begin at 0.2 Jybeam"^ and increase in steps of 0.1 Jybeam"^. 
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Fig. 9. — G31.41 channel maps of a) C^^O and b) H2S emission from 85 to 107 kms~^ increasing in steps of 2 kms~^. Contours begin at 
0.45/0.3 Jybeam"^ and increase in steps of 0.3/0.2 Jybeam"^ for C^^O and H2S respectively. Plus signs mark the location of water masers 
from Hofner & Churchwell (1996), the position of the UCHll region is marked by an open square and the filled circle marks the position of 
the 1.4-mm continuum peaJf. The beam is marked by the ellipse in the bottom left-hand corner of the image. 




Fic;. 10. G9.62 position-velocity diagrams along a) the core and b) the outflow axis respectively. The position angle of each cut is shown, 
measured positive east of north. The box in the top right-hand corner represents the effective resolution on each axis. 
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Fig. 11. — G10.47. position- velocity diagrams along a) the major and b) the minor axes respectively. The position angle of e3x:h cut 
shown, measured positive east of north. The box in the top right-hand corner represents the effective resolution on e3x:h axis. 
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Fig. 12. — G29.96. position-velocity diagrams along a) the core and b) the outflow axis respectively. The position angle of each cut 
shown, measured positive east of north. The box in the top right-hand corner represents the eifective resolution on each axis. 
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Fig. 13. — G31.41. position-velocity diagrams along a) the core and b) the outflow axis respectively. The position angle of each cut is 
shown, measured positive east of north. The box in the top right-hand corner represents the effective resolution on eax;h axis. 



